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Abstract 



Top quark mass suppressed terms are calculated for the virtual amplitude for Higgs 
production in gluon fusion at three- loop level, i.e. O(a^). The method of asymptotic 
)—} \ expansions in its automated form is used to evaluate the first three non-vanishing 

^ — , ■ orders in terms of Mg/M t 2 , where the first order corresponds to the known results of 

t-H \ the effective Lagrangian approach. 

Q-f 1 Introduction 

Radiative corrections to Higgs production through gluon fusion are known to be unusu- 
ally large [H [21 [31 0] . The inclusive next-to-next-to-leading order (NNLO) cross section 
, a(pp/pp — > H + X) exceeds the LO prediction by roughly a factor of two at LHC en- 

ergies, and even up to a factor of three at the Tevatron El [7]. Recent compilations 
of the currently available contributions to the production cross section can be found in 
Refs.pHIS]. 

The current NNLO prediction is based on the assumption that the top mass dependence is 
largely determined by the LO expression, while the higher order terms can be evaluated in 
Q^ \ the limit of infinitely heavy top mass Mt [W\ El HI [11]. At NLO, where a comparison with 

the full mass dependence of the cross section is possible, the heavy-top approximation 
is valid at the 2-3% level for Higgs masses Mh < 2M t (see, e.g., Ref. [E])- Even at 
Mjj ~ 1 TeV, the deviation from the full NLO result amounts to only about 10%. 

The fact that the heavy top limit works so well is at first sight surprising, because it 
assumes that M t is larger than any other scale in the process. This is certainly not 
the case at the LHC with a prospected hadronic center-of-mass energy of = 14 TeV. 
However, one can argue that since the cross section is dominated by soft gluon radiation 
parton scatterings with energies y/1 much larger than 2M t are strongly suppressed. 

It is indeed observed that an expansion of the partonic cross section a in powers of (1 — z), 
where z = Mu/s, converges rather quickly to the exact result [5]. On the other hand, 
resummation of the soft terms does not lead to a big effect at any of the three lowest 
orders in perturbation theory |13j . 

Recently it has been suggested that the size of the radiative corrections is due to the 
transition from space- to time- like momenta, and in fact, numerical studies show that the 
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bulk of the radiative corrections can be obtained by resumming the leading 7r 2 -terms that 
arise from this transition [14J. 

These unresolved issues leave one with a certain amount of doubt as to the use of the 
heavy-top limit at NNLO. There is however surprisingly little activity in the field that 
addresses the validity of this approximation. Besides the NLO calculations for the inclusive 
cross section mentioned before [lOj El [H [Tl] , there are studies concerning the mass effects 
on differential distributions [151 [To] [T7] which allow one to derive validity ranges on the 
kinematical variables. Furthermore, in Ref. |18j . the effects of the partonic high-energy 
region on the total cross section have been studied by deriving the leading behaviour in 
this limit. 

A rather direct way to check the heavy-top limit is to evaluate formally subleading terms. 
In this paper, we consider them for the purely virtual corrections at NNLO. While they 
do not correspond to a physical quantity, they constitute an important gauge-invariant 
ingredient to the full inclusive cross section. Note that at NLO, the virtual corrections are 
known in closed analytical form for arbitrary values of Mt \19\ [20l [2T] . 

Our approach is very similar to the calculation of the top mass suppressed terms to the 
Higgs decay rate into gluons, described in Ref. [22]. One might be tempted to use this 
result obtained for the decay rate as an estimate of the effects for the gluon fusion process. 
However, one should recall that the kinematics of the two processes are very different. 
In particular, the top quark mass is indeed the largest scale for the decay, so that the 
expansion in M^/Mt remains within the radius of convergence. This is not the case for 
the higher order corrections to the gluon fusion production process involving real radiation 
of gluons and quarks. The partonic center-of-mass energy \/i can well exceed the threshold 
value of 2M t , and a series expansion in the limit of large top mass becomes questionable. 
For the purely virtual effects though, which are the subject of this paper, the partonic 
center-of-mass energy is fixed to Mh which, according to the limits derived from electro- 
weak precision fits, can safely be assumed to be lighter than twice the top mass. They 
will therefore be a useful ingredient for any possible treatment of the full hadronic cross 
section, be it inclusive or exclusive. 

2 Method 

Sample diagrams that contribute to the virtual corrections to gluon fusion at LO, NLO, 
and NNLO are shown in Fig.[TJ An efficient and algorithmic procedure for evaluating them 
in terms of a consistent expansion in M^/M t is the well-known method of asymptotic 
expansions (see, e.g., Ref. [23]). In our case, it expresses the original diagrams as a sum 
of convolutions of massive vacuum with massless vertex integrals. The diagrammatic 
representation of this procedure is shown for two particular diagrams in Figs. [2] and El 

We generate the diagrams with the help of qgraf [24] and pass them to q2e/exp [25, 26j, 
which automatically carries out the expansion. The resulting 1-, 2-, and 3-loop vacuum 
integrals are evaluated by MAT AD [27]. For the 1- and 2-loop vertex integrals we use the 
method of Ref. [28] by applying the relevant modifications [29] to MINCER [30] . 

The colour and Lorentz structure of the physical amplitude is given by 

P£{qi,Q2) = $ ab (QvQ2 9^ - , (i) 
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(a) (b) (c) 




(d) (e) (f) 




(g) (h) (i) 

Figure 1: Examples of Feynman diagrams contributing to the virtual corrections 
at NNLO in the gluon fusion process. The solid lines denote top quarks, the 
springy lines are gluons, and the dashed line is the Higgs boson. In diagram (i), 
the bubble insertion can be a top quark or any other quark. 





Figure 2: Diagrammatic representation for the asymptotic expansion of a particu- 
lar Feynman diagram in the limit M| <C 4M t 2 . The dia grams left of ® represent 
subdiagrams of the original diagram that are to be expanded in the momenta 
corresponding to the dotted external lines before the loop integration. In this 
way, it is apparent that the original integral, depending on and M t 2 , is de- 
composed into products of "tadpole" integrals with vanishing external momenta 
and massless vertex integrals. The shaded blob in the diagrams right of (8) repre- 
sents an effective vertex given by the result of the diagram left of (8> (for details 
of asymptotic expansions, see Ref.[23], for example). The three terms right of 
"— are proportional to 7V t 3 , MfNu, and NtN%, respectively (cf. Eq. ^ below). 
Subdiagrams without external mass scales are not shown. 
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Figure 3: Diagrammatic representation for the asymptotic expansion of the Feyn- 
man diagram in Fig. [I] (i) in the limit Mg <C 4M t 2 , when the bubble insertion is 

a top quark. The second term on the r.h.s. is a source of the term ~ in 
Eq. ([9D below. 



where q^ and q% are the external gluon momenta, and a and b are the corresponding 

r>ab i 



colour indices. We contract the amplitude with P^(qi,q2) m order to arrive at a scalar 



expression in Lorentz and colour space. 

Before the massive two- and three- loop integrals are passed to MAT AD, we need to eliminate 
any external momenta in their numerators by appropriate decompositions into invariants, 
e.g. 



dl - = 5<W '») 7 d ' • • ■ -,/:■• „,»,- • (2) 

where the dots represent factors that are independent of I. 

There are two diagrams at one-loop level, 23 at two-loop level, and 657 at three-loop 
level, and the calculation of the l/M t 2 -suppressed terms takes about 5 • 10 4 s, with the 
computationally most expensive one shown in Fig. [31 



3 Results 

Before we present the results, let us introduce some useful notation. The renormalization 
scale n appears in our calculation only through the factors 

M h = ^ e {j£) 1 Mt= {wT M > ^ = ex P[ e (-7E + ln4^)], (3) 

with Euler's constant 7e ~ 0.577216 . These expressions are understood as their Laurent 
series in e = (4 — D)/2, where D is the number of space-time dimensions used in the 
calculation. 

The perturbative coefficients typically contain the transcendental numbers C, n = C( n )j 
where £ is Riemann's zeta function. The particular values occurring here are 

2 4 

C2 = ~ 1-64493 , Cs ~ 1.20206 , ^ = J- w 1.08232 . (4) 
o 90 
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Throughout this paper, bare quantities are labeled by a superscript "B" . Note that since 
the diagrams are evaluated with a spectrum of six quark flavours, renormalization has to 
be performed accordingly. To perform on-shell mass renormalization or conversion of a s 
from the six- to the five-flavour scheme one must keep the proper number of higher order 
terms in e due to the presence of infra-red poles. Furthermore, in order to arrive at a 
physical result, the external gluons must be renormalized on-shell. 

For convenience, the number of light flavours n\ is kept as a free parameter; the physical 
case corresponds to n\ = 5. The virtual cross section for the process gg — ► H can be 
written as 



7T 



1 
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5(l-z)\M t H(e)h(a»)\ 2 , 



(5) 



where 



H(e) = T(l + e) 



+ 



1 /M H Y 



3 1 , 
1 + ~e+ -e' 



168 \M*J V ' 2~ ' 2 
The amplitude is expanded in terms of a perturbative series: 

„B /^3 X 2 



7T y 7T J 



(6) 



(7) 



where the coefficients are functions of Mh, M t B , and the renormalization scale fi. In 
our approach, they take the form 



M H 

M t B 



M H 
M t B 



(8) 



The leading terms have been calculated in the framework of an effective Lagrangian. 
However, for consistency, we present them here in a form that is directly compatible with 
the mass suppressed terms to be presented below: 



h, 



4 2 ) =a (2) +c (2) +a (l) c (l)_ C (D,B a (l) ) 



(9) 



where 
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+ 0(e), 



5 



are the perturbative coefficients of the effective Higgs-gluon vertex as presented in Ref. [29J 
which we quote here for the sake of completeness. Furthermore, we find 

"3 
4 



„(2) 
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It may be worth noting that c^- 1 ' and correspond to the one- and two-loop results for 
the bare coefficient function of the effective LagrangianQ 
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Cf 



\— r(i + e) 
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a? id ( a® 
1 + — c\ ' + — 
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,(2) 



+ 0(ai 



(12) 



Finally, 



(13) 
6 to 



,,,+0(e 3 
6 \e 2" 3 ' 

is the 1-loop term of the bare decoupling constant for a s for the transition from n / 

rtj = 5 flavour QCD (see, e.g., Ref. [32]). The origin of the term involving Cg B i n Eq. Q 

is the fact that the coefficients oW in Ref. |29] were evaluated in 5-flavour QCD, while the 

of Eq. © are based on 6-flavour QCD. Therefore, diagrams like the right-most one in 

Fig. [3] do not have a correspondence in the effective theory calculation of Ref. [29|. 



The expressions presented so far correspond to known results and have been included in 
this paper only for the sake of the reader's convenience. They should facilitate any imple- 
mentation of the newly calculated terms to be presented below. Besides that, they serve as 
a useful check of our setup. It should be noted that in our approach, we directly calculate 
the coefficients hm\ and the decomposition into and is just for comparison to the 
literature. 



The new results of this paper are the contributions to the virtual 3-loop amplitude that 
are formally suppressed by powers of Mn/M t . In the notation of Eqs. ([7]) and dH), the first 
two subleading orders reaco 
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1 Renormalizing Cf according to Ref. [31] (where Cf is called C? ), one may derive the coefficient function 
Ci quoted in Eq. (3) of Ref. [H]. 

2 For the sake of brevity we insert SU(3) colour factors. The result for general colour factors can be 
obtained upon request from the authors. 
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where h>2 is the two-loop result for which the e 1 and e u terms can be compared to 



Ref. [10J, with full agreement, of course. Note that since the leading order amplitude 
has been factored out in terms of H(e), the NLO and NNLO expressions for h start only 
at 0(1/ e) and 0(1/ e 3 ) respectively. The leading poles are thus fully determined by the 
leading terms in 1/M t . 

There are a number of checks that we can perform on our result: (i) calculating the 
amplitude in an arbitrary covariant gauge, we find it to be independent of the gauge 
parameter; (ii) replacing the projector of Eq. ([T]) by 5 ' ffi^tfe,!/ leads to a vanishing result, 
which also checks gauge invariance; (iii) the poles of order a™ +1 /e k , k = 1, ... ,2n were 
checked against the general formula of Ref. [33J (with additional input from Ref. |34j for 
the a 3 /e terms; see also Ref. [35]) and found to be in full agreement. 

Another observation is that in the ratio of the amplitudes taken at time-like and space-like 
momenta, all mass effects cancel and the result is the one given by the general formula of 
Ref. [MIEZj. Specifically, 



(16) 



h(M%) 


2 


a(Ml) 






a(-Ml) 



where 



htf) = Z? s H(e)h(a 



(17) 



with the on-shell gluon renormalization constant Z® s [32j (see also Ref. [38] ) . The ratio on 
the left-hand side of Eq. (fTBT) can also be found in Ref. [29] . It is understood that the bare 



quantities M t B and af in Eq. (|16p are expressed in terms of their renormalized values [39] . 



4 Conclusions 



The top mass suppressed terms for the virtual corrections to Standard Model Higgs pro- 
duction in gluon fusion were presented through three-loop order. They are an ingredient 
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for the inclusive NNLO QCD cross section. We have subjected the result to various checks 
and found full confirmation. 

The next steps towards the top mass effects in the Higgs production cross section will be 
the evaluation of the mass suppressed terms in the real radiation amplitudes. We defer 
this problem to a forthcoming publication. 

Upon completion of this paper, we became aware of a similar calculation [3D]. We have 
compared our results and found full agreement. 
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